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We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk 23 1 . In total 25 lines are detected, including 
CO 7 = 5-4 through / = 13- 12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the 
CO rotational levels up to / = 8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution 
of the CO lines over the rotational ladder above 7 = 8 requires the presence of a separate source of excitation for the highest CO lines. We explore 
X-ray heating by the accreting supermassive black hole in Mrk 231 as a source of excitation for these lines, and find that it can reproduce the 
observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this 
model strongly overpredicts the hot dust mass in Mrk 23 1 . Our favoured model consists of a star forming disk of radius 560 pc, containing clumps 
of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to 7 = 8. X-rays from the accreting supermassive black 
hole in Mrk 231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in 
Mrk 231. The extraordinary luminosity of the OH* and H2O* lines reveals the signature of X-ray driven excitation and chemistry in this region. 
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1. Introduction 

Carbon monoxide (CO) is a fundamental tracer of interstellar 
molecular gas. However, since only the lowest 3 rotational tran- 
sitions are relatively easily accessible with ground-based tele- 
scopes, the diagnostic use of higher rotational levels is poorly 
developed. This hiatus in our knowledge is becoming acute now 
that high-7 CO observations of high-z galaxies are becoming 
possible. 

With this in mind, we have embarked on the Herschel 
Comprehensive (U)LIRG Emission Survey (HerCULES), an 
open time Key Project on the ESA Herschel Space Observatorjy 
(Pilbratt et al., this issue). The HerCULES project will estab- 
lish a comprehensive inventory of the gas cooling lines in a 
flux-limited sample of 29 (Ultra)luminous Infrared Galaxies or 
(U)LIRGs, using high spectral resolution observations with the 
Fourier Transform Spectrograph (FTS) of the SPIRE instrument 
(Griffin et al., this issue), combined with observations of the 
[Cn] 158 ^lm fine and the [Oi] 63 and 146 //m fines witfi PACS 
(Poglitsch et al., this issue). Key aims of HerCULES are the de- 
velopment of the diagnostic use of the gas cooling lines in local 
(U)LIRGs, and establishing a local benchmark for observations 
of high-z galaxies with the Atacama Large Millimeter Array. In 
addition, since tfie FTS yields full spectra, any otfier luminous 
emission lines detected (e.g., of H2O) will be available for study. 



* e-mail: pvdwerf@strw . leidenuniv . nl 

' Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



ISpaans & MefierinM(l2008h fiave sfiown tfiat X-ray excitation 
of tfie gas (e.g., by an AGN) and UV irradiation by young mas- 
sive stars produce very different luminosity distributions over tfie 
CO rotational lines. Pfiysically, tfie difference arises because X- 
rays penetrate a larger column density of gas tfian UV pfiotons, 
and are less effective in dissociating tfie molecules. In addition, 
wfiile tfie gas fieating efficiency in a pfioton dominated region 
(PDR) is less tfian 1 %, in X-ray dominated regions (XDRs) tfiis 
efficiency is 10-50%. As a result, for comparable irradiated en- 
ergies, XDRs tend to fiave larger column densities of warmer 
molecular gas tfian PDRs, and will produce mucfi more luminous 
emission in tfie fiigfi-7 CO lines. In contrast, PDRs are more ef- 
ficient tfian XDRs in fieating tfie dust. X-rays also give rise to 
significant ionization in tfie molecular gas and tfierefore drive an 
efficient ion-molecule cfiemistry, leading t o pronounced cfiem- 
ical d ifferences between PDRs and XDRs dMeiierink & Spaansl 
l2005h . Testing and using tfiese diagnostics is one of tfie principal 
aims of tfie HerCULES project. 

In tfiis Letter, we discuss tfie first results of tfie HerCULES 
programme and present tfie SPIRE FTS spectrum of tfie nearby 
ULIRG Mrk231 (UGC8058, IRAS F12540-H5708), tfie most 
luminous galaxy in tfie Revised IRAS Brigfit Galaxy Sample 
(Sanders et al. 2003). Adopting z = 0.042170 as tfie fiefiocen- 
tric redsfiift of Mrk 231, correcting for tfie local flow, and apply- 
ing a flat 5-year WMAP cosmology (//q - 70 km s"' Mpc"', 
Qa - 0.73) yields a luminosity distance Dl - 192 Mpc, 
witfi 1" corresponding to 0.856 kpc, as provided by NEEQ. 
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Fig. 1. SPIRE FTS spectrum of Mrk 23 1 . Line identifications are given in red for CO lines, in blue for H2O, in magenta for OH^ 
cyan for HaO^, and in green for the remaining lines. 
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The derived 8 - 1000 fim luminosity of Mrk 231 is then Lir - 
4.0 • 10'^ Lq. Mrk 231 contains a luminous, optically visible 
AGN, classified as a Seyfer t 1 or a Broad Absorption Line 
QSO (iBoksenberg et al.ll977h . A highly absorbed p ower-law X- 
ray spectrum was observed by iBraito et aP (l2004l) with Lx - 
6!° 3 ■ 10"*^ erg s-> between 2 and 10 keV. However, Mrk 231 
also contains a kpc size disk harbouring int ense star formatio n 
as shown by high resolution radio imaging dTavlor et al.lll999h . 
Interferometric imaging of CO 7 = 1-0 and 7 = 2-1 emis- 
sion shows an inner disk of radius ~ 520 pc, containing 45% of 
the total molecular gas mass, and embedded in a more extended 
and diffiise emission component; the total molecular gas mass is 



5 ■ lO'' Mo (iDownes & Solomonl[T99 8^. The lowest part of the 
CO lad der (up to / = 6-5), was analysed by Papadopoulo s et all 
(|2007|) . who showed that the integrated CO emission can provide 
a significant contribution to the total gas cool ing. Indications for 
X-ray-driven chemistry have been fou nd by Aalto et aP (l2007h 
in HNC and HCN line ratios, and by iGonzalez- Alfonso et al.l 
(2008|) in the abundances of OH and H2O observed in absorp- 
tion wit h ISO. In a comprehe nsive study of ULIRGs and low-z 
quasars. [Veilleux et al.l (l2009h derive a fraction of 70+15% for 
the AGN contribution to the far-infrared luminosity of Mrk 231, 
with the remainder coming from star formation. 
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2. Observations, data reduction and results 



Mrk 231 was observed in staring mode with the SPIRE 
FTS on December 9, 2009, as part of the Herschel Science 
Demonstration Program. The high spectral resolution mode was 
used, yielding a resolution of 0.04 cm"' over both observ- 
ing bands: the long wavelength band covering 14.9-33.0 cm"' 
{A = 671-303 yum, v = 467-989 GHz) and the short wave- 
length band covering 32.0 - 51.5 cm"' {A - 313-194 yum, 
V = 959-1544 GHz). In total 50 repetitions (100 FTS scans) 
were carried out, yielding an on-source integration time of 
6660 s. A reference measurement comprised of 120 repetitions 
was used to subtract the combined emission from the sky, the 
telescope and instrument itself. The data were processed and 
calibrated (using the asteroid Vesta) as described in Swinyard 
et al. (this issue). Since the CO extent of Mrk 231 is at most 2" 
dDownes & Solomonll 19981) . while the SPIRE beam varies from 
17" to 42" over our spectrum, calibration procedures appropri- 
ate for a pure point source were adopted, and no corrections 
for wavelength-dependent beam coupling factors were neces- 
sary. Because of the excellent match in the overlap region of the 
two spectrometer bands (32-33 cm"'), the bands were simply 
averaged in this region. 

The full SPIRE FTS specti'um of Mrk 23 1 is shown in Fig. [1] 
It shows a total of 25 well detected lines. The full CO ladder is 
detected with 9 lines from CO J - 5-4 to 7 = 13-12. In addi- 
tion, 7 rotational Unes of H2O are detected, and the two [C i] fine 
structure lines and the [N 11] fine structure line, as well as rota- 
tional transitions of CH^ and HF. Very surprising is the detection 
of luminous emission from OH^ and H20^ . While a possible de- 
tec tion of absorption in higher tr ansitions of OH^ was reported 
bv iGonzalez-Alfonso et al.l (12008.) . this is the first astronomical 
detection of H20^ except in comets. 

The lines are detected superposed on a continuum which 
rises towards the short wavelength side, and represents the 
Rayleigh-Jeans tail of the dust emission in Mrk 231. A small 
difference in the thermal background between the source and 
the reference observation manifests itself as a moderate flux ex- 
cess in the continuum at the longest wavelengths. Line fluxes 
are not affected by this artefact. For these early observations, no 
matching observations with the SPIRE photometer, to measure 
the continuum levels, were carried out. Pending this, quantita- 
tive analysis of the continuum is premature and in this Letter we 
restrict ourselves to the spectral lines. 

Before fitting line profiles, we subtracted the continuum 
emission using a grey-body fit made to the underlying spec- 
tral energy distribution (SED). Any remaining large-scale ripples 
were removed using a polynomial or sine wave fit. Line fluxes 
were recoved from this baseline-subtracted spectrum by itera- 
tively fitting model line profiles to this spectrum. These model 
line profiles are the convolution of the FTS full resolution instru- 
mental response (a sine function) with the underlying Gaussian 
line profile of the emission from the galaxy. The systematic 
uncertainty in the flux scale for the lines is 20-30% over the 
21-52 cm"' waveband, but significantly higher below 21 cm"' 
(which will improve when brighter calibration sources become 
available). We note that the RMS fluctuations in the spectrum 
are higher than the thermal noise, as a result of a fringe due to 
a standing wave in the instrument, which affects the accuracy of 
the derived parameters for the faintest lines in the spectrum. The 
removal of this fringe, together with a search for additional faint 
lines, is the subject of ongoing work. 
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Fig. 2. Luminosities of CO lines from Mrk 231. Filled symbols 
represent measurements from the SPIRE FTS spectrum, while 
ground-based measurements are denoted with open symbols. 
Coloured lines indicate two model PDR components (red and 
green lines) and an XDR component (blue line). The sum of 
these three components is indicated by the black Une and fits the 
CO measurements. In the legend, n denotes the number density 
of hydrogen nuclei (n = rtu + 2nu-,) in cm"^. Go denotes the in- 
cident UV flux in units of 1.6 • 10"^ erg s"' cm"^ for the PDRs, 
and Fx the incident X-ray flux for the XDR. The legend also 
indicates the relative emitting areas of the three components. 



3. Discussion 

3.1. CO excitation 

We combine the CO line fluxes from the spectrum shown in 
Fig. [U with ground-based measurements of the lower lines 
dPapadopoulos et aP l2007l and references therein) in order to 
construct the CO rotational excitation diagram shown in Fig.|2] 
It is seen that an approximately flat luminosity distribution is 
obtained for the lines from J - 5 upwards. Note that the CO 
J - 10-9 line is blended with the H2O 3 1,2-22,1 line which is 
expected to have some luminosity (1.8 ■ 10^ Lq for the model by 
Gonzalez- Alfonso et al., this issue), and this may account for its 
somewhat high flux. The total luminosity measured in the CO 
lines up to 7 = 13-12 is (2.8 + 1.1) • 10** Lq. Note that only 4% 
of the total CO line luminosity is contained in the lowest 3 tran- 
sitions. For compa rison, in our Milky Way this fraction is 43% 
( lFixsenetalJll999l) . 

The approximately flat distribution of CO line luminosity 
with rotational level indicates that several excitation components 
must be present, since individual components always produce a 
more peaked excitation diagram. We model th ese components 
using the one dimensional PDR/XDR models of lMeiierink et al.l 
(I2OO7 I). as shown in Fig.|2] The CO lines up to 7 = 8-7 can be 
produced by a combination o f 2 PDRs, in qualitative agreement 
with the decomposition by Papadopoul os et al](l2007h . However, 
a challenge is presented by the highest CO Unes: J - 13-12 and 
J = 12-11, arising from energy levels 503 and 461 K above 
the ground state. As shown in Fig. [31 these lines are strongly un- 
derproduced by the PDRs dominating the emission in the lower 
lines, since the resulting gas temperatures are not high enough 
for significant population of the 7 > 10 levels. These lines there- 
fore require the presence of a third excitation component, which 
can be either an XDR or a high excitation PDR. 
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A model fit with an XDR producing the highest CO lines is 
shown in Fig. |2] The required X-ray illumination for this XDR 
can be produced by the AGN in Mrk 23 1 ( Braito et al. 2004), out 
to a distance of 160 pc from the nucleus, ignoring absorption. 
The ratio of radiating surfaces in the model shown in Fig.|2]im- 
pUes an extended low excitation PDR component (green curve), 
with a less extended and denser central XDR region (blue curve). 
Dense clouds with a smaller surface, close to massive stars and 
probably embedded in the more diffuse component, account for 
the medium excitation component (red curve). 

Alternatively, a very dense, high illumination PDR can ac- 
count for the highest CO lines. A good fit is found with n - 
10^'^ cm""* and Go = lO'' and a surface ratio from medium to 
high excitation of 1.0 : 0.03. Here the small surface area for the 
high excitation PDR indicates a number of small high density 
clumps in a very strong UV field. Since the radiating surface of 
the high excitation PDR is about 30x smaller than that of the 
medium excitation PDR, but its density about 30x larger, the H2 
masses in these two components must be comparable. For an 05 
star, the required Go = lO*" is reached at a radius of 0.3 pc. For a 
star formation rate of 100 Mq yr"' and a power law initial mass 
function with slope -2.35 between masses of 0.3 and 120 M©, 
there are 7 .6 • 1 0^ stars of spectral type 05 or earlier in Mrk 231. 
The total volume with Go > 10^ is then 8.6-10'* pc^ , while (for a 
520 pc radius disk with 15 pc thickness, following ' Davies et al.l 
(120041) ) the total volume of the gas disk is 1.3 • 10^ pc^. In other 
words, in this scenario approximately half of the molecular mass 
would have to be contained in 0.7% of the total volume, and lo- 
cated within 0.3 pc from an 05 (or hotter) star. Efficient UV 
heating by the Go = 10^ radiation field would heat the dust in 
these clumps (in total ~ 50% of the total dust mass in Mrk 231) 
to a temperature of about 170 K. In contrast, in the XDR model, 
where dust heating wo uld be less efficient, the du st temperature 
would only be ~ 70 K (lMei)erink& SpaansllIOOSi) . 

These predictions can be tested by analysing the infrared 
SED of Mrk 231. Gonzalez- Alfonso et al. (this issue), found that 
the hot {T = 150-400 K) component in the SED of Mrk231 
accounts for about ~ 20% of the total infrared luminosity, but is 
produced by only 0.02% of the total dust mass. This result limits 
the fraction of gas within 0.3 pc from an 05 star in Mrk 231 to 
much less than the ~ 50% required to produce the highest CO 
lines with a high excitation PDR. This problem does not exist for 
the XDR model, which predicts most of the dust to be cooler 



the high excitation PDR, and the OH^ abundance is larger by a 
comparable amount. A similar argument can be made for H20^, 
which is formed by OH+ + H2 ^ H2O+ + H. 

The extraordinary luminosity of the OH^ and H20^ (and to 
a lesser extent CH^) lines in Mrk 231 is underlined by a com- 
parison with the SPIRE spectrum of the Orion bar PDR (Habart 
et al., this issue), which shows no trace of OH^ or H20^, and 
only weak CH^ emission, while in Mrk 23 1 the lines are only a 
factor 2-3 fainter than the CO lines. While enhanced cosmic ray 
fluxes in a starburst environment will increase the degree of ion- 
ization and hence the production of OH^ and H20^ in a PDR, 
they do not elevate the temperatures to the level required to pro- 
duce the highest CO lines ( Meiierink et al. 2006) . It is thus the 
combination of strong high-7 CO lines and high OH^ and H20^ 
abundances that reveals X-ray driven excitation and chemistry in 
Mrk 231. 



4. Outlook 

We have shown that the SPIRE spectrum of Mrk 231 reveals 
both PDR and XDR emission lines, and made a separation of 
these components. A key goal of the HerCULES project will be 
using this decomposition for a quantitative separation between 
star formation and black hole accretion as power sources for the 
infrared luminosities of dusty galaxies. In the case of Mrk 231, 
this issue will be addressed in a forthcoming paper (Meijerink 
et al., in prep.). Data will be obtained as part of the HerCULES 
program for an additional 28 objects, which will enable us to put 
the results presented here on a statistically significant footing. 
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3.2. XDR chemistry 

The extraordinarily luminous emission from the molecular ions 
H20^ and OH^ reveals the chemical signature of an XDR. 
Assuming that the emission arises from a disk with 160 pc ra- 
dius (as derived above), we can derive column densities in the 
upper levels of the relevant transitions, which results in values 
A^up ~ 3.0-3.5-10'^ cm"^ for both species. Modeling the nuclear 
molecular gas disk in Mrk 231 with a radius of 520 kpc and H2 
mass of 2.2 - 10^ M© then results in lower limits to the total OH^ 
and H20^ abundances relative to H2 of ~ 2 - 10"'° in the cen- 
tral 160 pc. Given the short radiative lifetimes (< 60 sec) of the 
upper levels involved, most H20^ and OH^ molecules will be in 
the ground state, and total abundances will exceed these lower 
limits by large factors. Such abundances require an efficient and 
penetrative source of ionization in the molecular gas, since the 
production of OH^ is mainly driven by H^ -H O — > O^ -H H fol- 
lowed by 0+ -H H2 ^ 0H+ + H and H+ -H OH -^ 0H+ + H. In 
the models considered in Sect. 13.11 the key species H^ is 2-3 
orders of magnitude more abundant in the XDR model than in 
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